Vortex structures of a two-component BEC for large anisotropics 
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We calculate the vortex structures of an elongated two-component Bose-Einstein condensate. 
We study how these structures depend on the intra-component and inter-component interaction 
strengths. We present analytic and numeric results respectively at weak and strong interactions; 
finding lattices with different interlocking geometries: triangular, square, rectangular and double- 
core. 
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I. INTRODUCTION 

One of the most exciting recent developments in cold 
atom experiments has been the production of artificial 
gauge fields, which couple to neutral atoms in the same 
way that magnetic fields couple to charged particles [1, 2]. 
While the greatest excitement surrounds the possibility 
of producing analogs of fractional quantum Hall states 
or topological insulators [3], these experiments also al- 
low one to perform some extremely interesting experi- 
ments on vortices in Bose-Einstein condensates. These 
structures are particularly rich for multi-component gases 
[4-7]. A key difference between the latest NIST exper- 
iments, and previous ones on rotating condensates [8- 
13], is that the techniques used to generate the artificial 
gauge fields naturally lead one to consider very narrow 
geometries. Here we theoretically investigate the vortex 
structures in a two-component Bose-Einstein condensate 
confined to a relatively narrow channel. We find a rich 
phase diagram which complements our understanding of 
the isotropic 2D system [4-7]. 

In the NIST experiments the hyperfine states of 87 Rb 
are coupled by a series of Raman lasers. In the dressed 
state picture [1, 2], taking the strong Raman field limit, 
one arrives at an effective single particle Hamiltonian 

H, = [ Px - A x {y)] 2 /2m+pl/2m + p 2 z /2m + V cS . (1) 

The effective gauge field is bounded, —hkji < A x {y) < 
hkn, where hkn is the momentum kick from absorbing a 
photon from one beam and releasing it into another. One 
can generate arbitrarily large magnetic fields B = d y A Xl 
but due to the bounded nature of A x , these can only ex- 
ist over a finite extent in the y direction. It is therefore 
natural to choose the effective potential V e fi(y, z) to re- 
strict the particles to a small region of y, z, modeling it 
as V e ff(y, z) = muj* 2 y 2 /2 + muj z 2 z 2 /2. V c s includes both 
the external potential and one induced by the Raman 
lasers, and w* is an effective trap frequency. 

The impact of the magnetic field on a single compo- 
nent Bose-Einstein condensate in such a narrow channel 
has been extensively studied [17-19]. As one increases 
the magnetic field from zero, the Bogoliubov spectrum 



starts to develop a "roton" minimum, whose energy falls 
with increasing B. When the roton energy hits zero, the 
system becomes unstable to undulations. At higher mag- 
netic fields vortices enter the channel. Here we extend 
these results to the two-component gas. 

The vortex structures in isotropic 2D two-component 
condensates are quite rich [4, 5]. Depending on the 
relative strength of the inter-component and intra- 
component scattering one can find interpenetrating lat- 
tices with different geometries: triangular, square, rect- 
angular and double-core. Additionally, one sees inter- 
esting hysteresis [5]. Wc find similar results in the 
anisotropic geometry. 

A final motivation for thinking about the role of mag- 
netic fields in very narrow geometries comes from solid 
state physics. In 2005, Seidel et al [16] showed that the 
quantum Hall effect in narrow tori can be connected to 
charge density waves, with the fractionally charged vor- 
tices mapping onto fractionally charged domain walls. In 
this paper we will focus on the Bose condensed regime, 
and will not be able to comment on this interesting 
physics. We do note, however, that our system undergoes 
a charge density wave instability before vortices enter the 
system. Investigating this instability in the low density 
"quantum" limit may be fruitful, even though the physics 
is likely to be very different from that in [16]. 



II. MODEL 

There are four important lengths in Eq. (1), d z = 

\Jh/mui z , d y = ^Jh/mui*, £ = \J h/B, and R = hkjj/B, 

where B is the peak value of B. These correspond to the 
trap length in the z and y direction, the magnetic length, 
and the spatial range over which the effective magnetic 
is nonzero. If d z is sufficiently small, the extent in z di- 
rection is sufficiently compressed and the 3D gases can 
be treated as 2D cloud. As we argue below, a sufficient 
condition will be d z < £. If d y <C R, the cloud will be 
confined to a region where the vector potential varies lin- 
early and we can approximate: A x = By. In terms of 
the raw experimental parameters, B — 5'hkji/D,, where 
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f2 is the Raman Rabi frequency, and 6' is the gradient of 
two-photon detuning along y direction. 

Taking *i(r) and *2(r) to annihilate atoms in the two 
(pseudo)-spin states, the short range interactions will be 



H mt = \J <M.9l*l*l*l*l +.92*2*2*2*2 



+2^12*1*^*2*1 



(2) 



The coupling constants are related to scattering lengths 
a s via g — Anh 2 a s /ra. We will consider the case g\ — 
.92 = .9, and gn = g m - For most experiments g\ w g-i ~ 
512, but one gains insight by relaxing this condition. 

Following Sinha and Shlyapnikov [17], we diagonalizc 
the single particle Hamiltonian in Eq. (1). The eigen- 
states are labeled by three quantum numbers K,n,n', 
with energies 



E nn ' (K) = £K 2 + nhuj z + n'M) c 



(3) 



B/m, and £ = 



where uj 2 = uj* 2 + uj 2 , uj c 
H 2 u}* 2 /imuj 2 £ 2 , I = \Jhjmujc, and we have neglected the 
zero-point energy. The continuous variable K = y/2£k 
is a dimensionless label of the momentum k along the x 
direction, while n and n' are discrete quantum numbers 
corresponding to the number of nodes in the z and y di- 
rections. If the interaction energy per particle (H int /N) 
and the characteristic "kinetic energy" (£K 2 ) are small 
compared to huj c and huj z one can truncate to the single 
particle eigenstates with n = n' = 0, which are of the 
form 



<f> K (x,y,z) = 



1 



(v-vk) 



(4) 



where yx = V2lu c K£/2uj c and L is the length of x direc- 
tion. 

By taking the system sufficiently dilute, it is easy to 
arrange (Hi nt /N) <C Huj c ~ fnv z . The other condition, 
£K 2 <C hu>c, will be valid at strong magnetic fields. 
For our vortex lattices, we find the characteristic di- 
mensionless wave- number to be K ~ 1, thus requiring 
uj* <C lo c . Combining this with the previous constraint, 
uj* ^> muj^/hk 2 ^, we see that our approximations break 
down unless the magnetic length is much larger than the 
wavelength of the Raman lasers, ie. uj c <C Kk 2 R jm. This 
also establishes our requirement d z < £. 

Letting ax annihilate the state in Eq. (4), the iV-body 
Hamiltonian is 

H/e = J2 r2 [ a U a ^ + 4^] + 1: E F t + («) F t(9) 

K q 



where 



S q - Kl -K 2 e I<2)2 a Kia a K2T . (6) 



if i if 2 



We take the continuum limit ^2 K —> {^/2L/Att£) J dK, 
and 5k — > (2^/2tt£/ L)5(K). The effective ID parameters 
are f3 = Nmguj 2 £/TrLd z uj* 2 h 2 , (3 m = (3g m /g. From these 
definitions of the effective interaction parameters /3,/3 m , 
we sec that increasing g, g m has the same effect as increas- 
ing the magnetic field B, increasing the confinement w z , 
or reducing the confinement uj*. 

We consider a variational wavefunction corresponding 
to a Bose-Einstein condensate which is periodic in the x 
direction, with dimensionless wavelength 2tt/Kq, corre- 
sponding to a physical wavelength A = 2\[2-k£I Kq 



exp 



cj \ n 



(7) 

Where a =t,| for two-component BEC. N a is the num- 
ber of particles in state a, and N = iV-j- + N±. The coef- 
ficients C na are normalized to Y^ n \Cna\ 2 = 1. In place 
of such a coherent state ansatz, some authors prefer to 
work with a "Fock" state 



H')f = 



1 



of 



E c «i«UJ i vac >- ( 8 ) 



For the quantities we are interested in, and \iP)f are 
equivalent, and the variational parameters C na have the 
same meaning in each case: |C ntT | 2 is the fraction of spin 
a particles with momentum uKq. 

We take = N±, and further restrict ourselves to 
considering symmetric or antisymmetric wavefunctions: 
C n ,o = iC-n.cr. We classify our state by the value 
of Kq and £ the number of non-zero C„'s which are 
needed to minimize the energy. For example, as illus- 
trated in Fig. 1, for (/3,/3 m ) ~ (3,1), we need only one 
n for each component: n = 0. We refer to this state 
as (£ti£j.) = (!)!)• This should be contrasted with the 
case at (/3, f3 m ) ~ (5, 2.5), where the energy is minimized 
by taking n = —1,0, 1 in both components, a state we 
label as = (3,3). Analytic expressions for the 

energies with small numbers of components are given in 
appendix A. 



III. RESULTS AT SMALL /3, f3 m 

The number of expansion parameters, (£-(■, grow 
with the magnitude of f3,0 m . For small j3 and j3 m , we 
only need a small number of terms in our wavefunctions, 
and we may use the analytic expressions in Eq. A1-A15 
to find the lowest energy state. We obtain a series of 
phase, as demonstrated in Fig. 1(a). 

At /3 rn = 0, the f and | atoms decouple, and the 
physics is identical to the single component case, £f = 
£l = £■ When j3 < 9.8, the ground state wavefunction 
has only one term (£ = 1), which displays a Gaussian 
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FIG. 1: Top panel: phase diagram at small dimensionless 
intra-species and inter-species interactions P,Pm- Larger 
P, p m corresponds to larger g,g m , larger magnetic field, larger 
confinement u z , or weaker confinement lo*. Hatched patterns 
represent states described by different number of Fourier com- 
ponents in each spin state: (ff, £4.). Bottom panel: the density 
profiles of the two-component wavefunction in the correspond- 
ing regimes. A color key for the density patterns is shown in 
Fig.2. 



shape along y direction. As j3 increases past j3 = 9.8, the 
£ = 1 wavefunction becomes unstable, and it undergoes 
a second-order phase transition to a density wave with 
£ = 3. As /3 increases to /J = 10.8, a first-order phase 
transition occurs to a state with £ = 2, characterized by 
a single row of vortices. These results have been exten- 
sively studied in Ref. [19]. 

In the regime of repulsive inter-component interac- 
tion (f3 m > 0), the j" and J, particles try to avoid each- 
other. For strong repulsive interaction {(3 m > /3), the two 
components undergoes a microscale phase separate [14], 
which needs large (presumably infinite) xi to describe. 
For weak repulsive interaction (J3 m < ft) rich structures, 
illustrated in Fig. 1 develop. For some (3,(3 m , we find 
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FIG. 2: The density profiles of two-component BEC at large 
P,Pm, where x,y are in the unit of y/2£. For (a) — (e), the 
parameters are p m /P = 0.1,0.5,0.8,1,1.5 respectively, with 
P = 1000. A color key is shown at the top. 



£t £l- For example, at (/3,/3 m ) = (8,2), £ T = 1 and 
£l = 2, and the 4- atoms have a row of vortices, while the 
f atoms show no structure. Under these circumstances 
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there is a degenerate state with the t an d I wavefunc- 
tions reversed. For other ft,ft m , there is a symmetry be- 
tween the two components with ^ = The vortices 
or density corrugations are displaced by half a period 
so that the density maxima of the f atoms line up with 
the density minima of the others. For example, when 
(ft,ft m ) = (10,3), each component displays a single row 
vortices, and the wavefunctions arc related by a transla- 
tion. 

In the regime of attractive inter-component interaction 
(ft m < 0), one wants to maximize the overlap of the j" and 
1 wavefunctions. Generically, this means that the wave- 
function of each component is identical, and the problem 
reduces to the single component case, but with a renor- 
malized interaction ft — > ft + ft m . The phase diagram 
for ft m < can be calculated from the phase diagram at 
ft m = by mapping each point: (ft,ft m ) -> (ft - /3 m ,0). 
For strong attractive interaction (—/3 m > ft), the BEC is 
unstable and expected to collapse, similar to the case of 
a single component BEC with attractive interaction [15]. 

IV. RESULTS AT LARGE ft, ft m 

For large ft and ft m , the analytic expressions for the 
energy become unwieldy. We numerically minimize the 
expectation (H), varying {C ncri K Q } in Eq. (7). The en- 
ergy landscape has many local minima, and we use a 
range of starting parameters to try to find the absolute 
minimum. We cannot rule out the existence of even lower 
energy states. Moreover, for some parameters we found 
that the energy differences between competing minima 
became extremely small. In an experiment it is doubt- 
ful that one would find the true minimum energy state. 
Rather than systematically exploring the large ft physics, 
we simply show a few examples. 

Our results, illustrated in Fig. 2, are richer than those 
seen in the single component gas [19], showing struc- 
tures similar to those in isotropic 2D studies [4, 5]. 
When ft rn <C ft one finds two interlocking triangular 
lattices, as in Fig. 2(a) where ft m /ft = 0.1. As one 



increases ft m the lattice structure changes: interlock- 
ing square/rectangular lattices are shown in Fig. 2(b,c) 
where ft m /ft = 0.5, 0.8. When ft m > ft the non-rotating 
system would be expected to phase separate. Here, one 
finds more intricate vortex structures at ft m ~ ft. At 
ftm/ft = 1, we find double-core vortices (cf. [5]), as in 
Fig. 2(d). At ftm/ft = 1-5 we find stripes, which are a 
microscopic version of phase separation. Close inspection 
of the image in Fig. 2(e), shows vortex cores in the low 
density regions. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the vortex structures in a two- 
component BEC with an artificial magnetic field, in an 
elongated geometry. Compared to the single component 
gas, the two-component vortex structures are more intri- 
cate. 

To experimentally investigate these structures, one 
needs to find a system where the interspecies interactions 
can be tuned relative to the intraspecies. One promis- 
ing approach is to use different atomic species for the 
two (pseudo)-spin states, and take advantage of an in- 
terspecies Feshbach resonance [20]. Some of the other 
newly condensed atomic systems may also be favorable 
[21]. If one cannot separately tune ft and ft m , one can 
still change the magnitude of them, fixing ft m / ft. As seen 
in Fig. 1 , such a cut through the phase diagram can still 
be quite rich, especially if ft m / ft ~ 0.25. 
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Appendix A: Analytic results for small number of 
components. 

Here we give analytic results for the energies of the 
states defined in Eq. (7), truncating the n-sums for each 
cr, and assuming symmetry/antisymmetry about the ori- 
gin for each component. 

State (£-(-, £4,) = (1, 1) is unique, and has energy 



For state (£4-, £4,) = (2,2), there are two distinct ex- 
tremal states. Both components could be symmetric 
about the origin (Ci-f = C_i t = C14. = CI 14 = l/y/2), 
or one of them could be antisymmetric (Ci-j- = C_i-|- = 
Ci_i = —C-i_i = 1/V2). This gives energies 



^2,2 



K% + \{B + B m ) + ^{8 + B m )e~ K " (A8) 



ez 



-8 e~ K o 



(A9) 



For state (£4-, £4,) = (2,3), the f component can be 
symmetric or antisymmetric, resulting in energies 



Kl + 8(2e~* K ° - e~2 A 'o - l) 
i/J m (l_2e-**3+2e-*oy 



-- K 

1 "n 



(A10) 



[t - Ae -kKl +2e -hKl +e -Klyi 



^2,3 



E. 



2,3 



8 m e 



(AH) 



For state (£4-, £4,) = (3, 3), the optimal wavefunction is 



£1,1 = \{8 + B m ) 
State (f t> £,) = (1,2), with 



(Al) 



V-t = y/1 ~ 2\e\ 2 4> Qt + ie(0jr ot + 0_jf ot ) (A12) 
^ = Vl-2|e| 2 0o; ± ie(4>Koi + 4>-k i) (A13) 



ipt = <Pot 



i»X J. + <P-K i) 



is unique up to translation, and has energy 

^, 2 = ^ 2 + ^(e-< + 5) + 5/3 m e-^o 

State (£ t ,£ ; ) = (1,3), with 
ifo = <fof 

V-i = VI - 2|e|2^ - e {4> Koi + <j>- Ko ±) (A6) 
has energy 



(A2) 
(A3) 

(A4) 
(A5) 



with e real. 

If 8 > 8 m > 0, the negative sign has lower energy, 



^3 = + 



2Kl + 8{-2e"^ K " +4e 



rA',1 



2) 



+ /3 m (2e-^_2)e 2 



+ 8(2e- K » + 4e 



8e-3 K " +3)e 4 (A14) 



+ B m {2e- K ° -4e-i R o +3)e 4 



If 8 > 0, 8 m < 0, the positive sign has lower energy, 



E lt3 = -(8 + p m ) 



Kq + 8{2e-^ K « - er* K « - 1) + 8 m {eT* A ° - 1) 



+8 (e~ K o + 2e~ 1 * K o - Ae~i K 'o + | ) 



(A7) 



E 33 = E 3 3 + 4,B m {e- K o/ 4 - e- K "' 2 ){e 2 + 2e 4 ) (A15) 



